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The performance of a scheme proposed for automated routine monitoring of 
deep-space missions is presented. The scheme uses four different tones (sinusoids) 
transmitted from the spacecraft (S/C) to a ground station with the positive iden- 
tification of each of them used to indicate different states of the S/C. Performance 
is measured in terms of detection probability versus false alarm probability with 
detection signal-to-noise ratio as a parameter. The cases where the phase of the 
received tone is unknown and where both the phase and frequency of the received 
tone are unknown are treated separately. The decision rules proposed for detect- 
ing the tones are formulated from average-likelihood ratio and maximum-likelihood 
ratio tests, the former resulting in optimum receiver structures. 


I. Introduction 

It has been proposed that automated routine monitoring of deep-space missions be provided by trans- 
mitting one out of n (typically n = 4) different subcarriers (tones) from the spacecraft (S/C) and then 
using a small automated terminal (for example, a 6-m low Earth orbiter terminal (LEO-T)-class) ground 
station to detect the presence or absence of each possible tone. The positive identification of each of 
the tones at the receiver will indicate different stages of the S/C, for example, S/C healthy, S/C needs 
help, S/C is going to transmit telemetry, etc. Since each of these tones is transmitted from the S/C to 
the ground over an additive white Gaussian noise (AWGN) channel along with the added possibility of 
Doppler distortion, the above-mentioned detection problem to be solved at the receiver can be formulated 
as a binary hypotheses test of signal plus noise versus noise only. In the most general ease, the signal 
is modeled as a constant power sinusoid with unknown [i.e., uniformly distributed on (— 7r,7r)] phase 
and unknown (i.e., uniformly distributed in some interval {fi,//) governed by the amount of Doppler) 
frequency. 

The optimum solution to problems of this nature is based upon maximum-likelihood (ML) consid- 
erations of the type discussed in VanTrees [1]. In particular, the solution takes the form of a binary 
hypothesis likelihood ratio test against a threshold whose value depends on the specified false alarm and 
detection probabilities, the available signal power- to- noise spectral density ratio, and the duration of 


1 This work was performed under a NASA Summer Faculty Fellowship at the Jet Propulsion Laboratory, Communications 
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the observation of the hypotheses. We shall see that there are, in principle, two detection/estimation 
philosophies suggested by the ML approach, corresponding respectively to what is commonly known as 
noncoherent detection, wherein no attempt is made to estimate the unknown parameters prior to detec- 
tion, and pseudocoherent detection, wherein an attempt is made to first estimate the parameters (using 
an ML approach) and then to use these estimates to aid in the detection process [2]. Since there appears 
to be some question about which is the better approach, we shall consider both approaches, discuss their 
philosophical differences, and compare their performances. 

This article is organized in two parts. In Part 1, we consider the problem of optimally detecting 
a sinusoidal signal of known amplitude (power) and frequency but of unknown phase [i.e. , uniformly 
distributed on (-7r,7r)] transmitted from a S/C to the ground over an AWGN channel. In so far as 
the optimum receiver design is concerned, the problem will be formulated as a binary hypothesis test of 
signal plus noise versus noise only with a single unknown parameter (i.e., carrier phase). In Part 2, we 
consider the added possibility of Doppler distortion, which produces an uncertainty in the received carrier 
frequency. Once again, the problem can be formulated as a binary hypothesis test of signal plus noise 
versus noise only, where now the signal is modeled as a constant power sinusoid with unknown phase and 
unknown frequency. Unfortunately, however, the theory for this case is not a s well developed in [1] as for 
the case where frequency is known. Nevertheless, other researchers [3-6] have examined problems of this 
type in the context of frequency-hopped (FH) or direct sequence (DS) spread spectrum communication 
systems, and we shall make use of their results wherever appropriate. 


Part 1. Known Frequency and Unknown Phase 
II. The Average-Likelihood Ratio Test 

A. Derivation of the Optimum Decision Rule and Associated Receiver Structure 

Consider the transmission of a fixed (known) amplitude sinusoid with known frequency and unknown 
phase over an AWGN channel. As such, the received signal can be modeled over the interval of observation 
0 < t < T as corresponding to either of two hypotheses, namely, 


r(t) — s(t, 9) + n(t) = V2P cos(uj c t + 0) + n(t) 


(la) 


when indeed the signal was sent (hypothesis H 1 ) or 


r(t) = n(t) 


(lb) 


when the signal was not sent (hypothesis H 0 ). In Eq. (la), P,u) c respectively denote the known signal 
power and radian carrier frequency, and 6 denotes the unknown carrier phase assumed to be uniformly 
distributed in the interval (— 7 r, 7 r). Also, n(t) denotes the AWGN with single-sided power spectral density 
No W/Hz. 

The optimum detection of a signal transmitted over an AWGN channel is the solution to the problem 
of finding the likelihood ratio (LR), defined as the ratio of the conditional probability density functions 
(pdf’s) of the received signal under the two hypotheses, namely, 


A (r(t)) - 


P(r(QlJTi) 

p(r(t)\H 0 ) 


( 2 ) 


and then comparing this ratio to a suitable chosen threshold to decide between H 1 and Ho , i.e., 
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( 3 ) 


Hi 

A (r(t)) > V 

H o 


In the case where all parameters of the signal are known, the evaluation of the numerator and denominator 
of Eq. (2) is straightforward, namely, 


p(r(£)|tfi) 


exp 


_1_ 

Wo 


T 


J (r(t) - s(t)) 2 dt 
0 


p(r(*)|tf 0 ) 



( 4 ) 


When the signal has an unknown parameter, e.g., the phase 0 , then to compute the numerator of Eq. (3), 
we must first condition the pdf p(r(t)\Hi) on the unknown parameter (9) and then average over this 
parameter, i.e., 


P 


(r(t) \Hi) = f 

J —IT 


p(r{t)\Hi,6)p e {6)d0 


( 5 ) 


where po(0) denotes the pdf of the unknown parameter 0. In our situation, the phase is assumed to be 
completely unknown and, thus, po(0) is a uniform distribution. Also note that this conditioning on the 
unknown parameter is now necessary in the denominator of Eq. (3) since the signal does not explicitly 
appear in p(r(t)\Ho) [see Eq. (4)]. Hence, combining Eqs. (3) through (5), the average-likelihood ratio 
(ALR) 2 becomes 


A(r(*)) - 


^ P(r(tm , D)dO S S’, “P {- ^ So 'MO - s(t, de 

p(r(tVH„) 1 f 1 r r 


J 

r pt) 

{ 1 . 

r j 

f 2 f T ( 

= exp <| 

{ ^oJ 

[ 2tt J 

1 exp < 

{No Jo r 

J 

r pt] 

{ 1 

r j 

{ 2 V 2 F r 

= exp j 

{ N 0 ] 

1 27 rj 

' exp < 

0 

£ 

✓ 


VtsN o 
r(t)s(t,9)dt ^ d6 


i 


(6) 


2 We shall refer to this formulation as an average-likelihood ratio (ALR) test to distinguish it from another (in general, less 
optimum) approach to be discussed shortly, in which a best (maximum-likelihood) estimate of the unknown parameter is 
obtained first and then used in the detection process. We shall refer to the latter approach as a maximum - likelihood ratio 
(MLR) test. This vernacular is not standard in the literature. What is important to understand here is that the words 
average and maximum refer to the manner in which the unknown parameter is handled, i.e., the estimation part of the 
problem and not the manner in which the decision on the hypothesis is made, i.e., the detection part of the problem. We 
shall be more explicit and mathematically precise about these differences later on in the article. 
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In arriving at the final result in Eq. (6), we have noted that the term exp j — / Q r r 2 (t)dt/No j is common 
to both the numerator and denominator and, thus, cancels, and also that 


exp \-'ki s2mdt }" exp {~^i c«v t < + s)4=»p{-f} 


( 7 ) 


assuming uj c T >> 1, as is typically the case. Defining the in-phase (I) and quadrature (Q) correlations 


[ r(t)y / 2 
Jo 


cos uj c tdt 


f r(t)V 2 
Jo 


sin uirtdt 


then Eq. (6) can be rewritten as 


A W f» = «P } £ /__ «P { ^ i “»(« + M } <» = <*P f - ^ } /o t 


(8) 


where 


l ± v'lfTif 


A _ 1 L s 
a = tan — 

-C«c 


( 9 ) 


Comparing A(r(t)) to a threshold rj is equivalent to comparing In A(r(t)) to In?;. Thus, taking the natural 
logarithm of Eq. (8), we obtain the equivalent decision rule 


Ho 


( 10 ) 


Finally, since In J 0 (x) is a monotonic function of its argument, x, and since PT/Nq can be absorbed into 
the decision threshold, then the decision rule of Eq. (10) can be further simplified to 


H x 

^1 L > < 
N 0 L < * 

Ho 


( 11 ) 


or, equivalently, 
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( 12 ) 


H i 

r2 > ,2^0 A 

< ? 4P 7 

H 0 


i.e., the optimum decision of signal present versus signal absent is determined from a comparison of 
the output of a square-law envelope detector with a normalized threshold, 7, whose value is determined 
from the specifications on false alarm probability and detection probability (see the next section). An 
implementation of the decision rule in Eq. (12) is illustrated in Fig. 1. 



Fig. 1. Average-likelihood (noncoherent) detector for detection of a single sinusoidal 
tone with known frequency and unknown phase in AWGN. 


B. Performance (Receiver Operating Characteristic) 

The performance of the receiver in Fig. 1 is described in terms of its false alarm probability ( Pf ), 
defined as the probability of deciding Hi (signal is present) when indeed Ho is true (signal is absent), and 
its probability of detection (Pd), defined as the probability of deciding Ho (signal is absent) when indeed 
Hi is true (signal is present). These probabilities are readily computed from knowledge of the first and 
second moments of the Gaussian random variables L c and L s [see Eq. (8)] under the two hypotheses, 
namely, 


Ho : E{L C } = E{L S } = 0 


var {L c } = var {L s } 


N 0 T 


H x : E{L c \e} = JPT cos 9 

E{L s \e } = -VPT sin 6> 


var {L c } = var {L s } = 


JVnT 


(13) 


To compute P F , we observe that, under hypotheses H 0 , L is a Rayleigh random variable (L 2 is a central 
chi-squared random variable). Thus, 
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* - pr[H ' w - priL 2 > - r /" (- 4 =) ^ 

= /j^ 2fiexp( - R2, " = “ p (-^) <i4) 

Similarly, we observe that, under hypothesis # 1 , L is a Rician random variable (L 2 is a noncentral 
chi-squared random variable). Thus, 


P D - - 1M* > •**.> - f (-£^1) /. (^) it 


/? 2 = (£{L C |0}) J + (£{L S |0}) 2 = PT 2 


> (15) 


J \j2~r/N,iT \ Z 


I 0 (Rd)dR = Q ( d, 


27 

AT n T 


where 


2 a 2PT = 2 E 
No ~ No 

is the detection signal-to-noise ratio (SNR) and Q(a,) 3) is the Marcum Q-function defined by [1]: 

t°° ( z 1 + q 2 \ 

Q(a,P)~ I 2 exp f j I 0 (az)dz 


(16) 


(17) 


Combining Eqs. (14) and (15) and eliminating the normalized detection threshold, one obtains the receiver 
operating characteristic (ROC) given by 


P D = Q (d, y/—2 In Pp) (18) 

which is illustrated in Fig. 2 for several values of the parameter d (or, equivalently, E/No). Alternatively, 
Pd I s plotted versus d 2 with Pp as a parameter in Fig. 3. 


III. The Maximum-Likelihood Ratio Test 

A. Derivation of the Optimum Decision Rule and Associated Receiver Structure 

Although the exact evaluation of the numerator of the likelihood ratio in Eq. (2), i.c., p(r(t)\H\) is 
obtained from the law of conditional probability as described by Eq. (5), namely, conditioning on the 
unknown parameter and averaging its distribution, it is also possible to approximate this numerator 
by first finding the ML estimate of the unknown parameter and then substituting this value into the 
conditional probability p(r(t)\H\,9). That is, we approximate p(r(t)\Hi) by p ^r(£)|//i , * in which 

case the likelihood ratio test (now referred to as the maximum-likelihood ratio (MLR) test) becomes 
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Fig. 2. ROC: frequency known and phase unknown. 


A (r(t)) ~ 


H i 

p{r{t)\H u 9 ML ) > 
p(r{t)\H 0 ) < 


Ho 


(19) 


We refer to this approach of first optimally estimating the phase and then using this estimate to aid 
the detection process as pseudocoherent detection. It is important at this point to emphasize that in the 
general context of problems of this type, i.e., detection of signals with completely unknown parameters, 
the performance of a receiver derived from MLR considerations (e.g., a pseudocoherent receiver) is never 
better than the performance of the ALR receiver (e.g., a noncoherent receiver ), which is indeed optimum 
under the assumed conditions. Thus, at best , one could hope that the MLR receiver would perform equally 
well as does the ALR receiver. In the next section, we shall indeed reveal the extent to which this equality 
in performance can be achieved for the problem at hand. First, however, let us derive the ML estimate 
of phase, namely, 9 ml, to be used in approximating the numerator of the likelihood ratio. 


The ML estimate of 6 is defined as 


9 m r = max 


p{r(t)\Hi,6) 
V p(r(t)\H 0 ) 


( 20 ) 


Using Eq. (4) in Eq. (20), it is straightforward to show that 


0 ml = max 
e 


exp {^oi> ^m*’*)*! = 


r 2V2 p p 

max exp < — - — / r(t) cos (a ) c t + 0) 
Q [ J o 


dt 


= max exp < ^^-Lcos(0 -F a) 
e iV 0 


( 21 ) 
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Fig. 3. Detection probability {Pq) versus detection SNR (d 2 ): 
(a) frequency known and phase unknown and (b) frequency 
known and phase unknown (expanded view). 


where the envelope, L, and the phase, a, are defined by Eq. (9) together with Eq. (8). Since L is 
positive and independent of 0, then the maximization required in Eq. (21) is achieved when the argument 
of tiie cosine function is equal to zero. Thus, 


L = 


( 22 ) 


An implementation of this ML estimator of the unknown channel phase is illustrated in Fig. 4. Also 
illustrated in Fig. 4 is the pseudocoherent detector that employs this ML phase estimator, which can be 
obtained by taking the natural logarithm of Eq. (19). We now find the decision rule based on the MLR 
test in Eq. (19) and compare it with that of the previously discussed ALR test. Using Eq. (22) in Eq. (19) 
gives, by analogy with Eq. (8), 




MAXIMUM-LIKELIHOOD PHASE ESTIMATOR 


Fig. 4. Maximum-likelihood phase estimator and pseudocoherent detector. 


r pt) f 2\/P /- ( PT) f 2\/P 

A(r(*)) exp exp Lcos ( 0 ml + «J | = exp / exp 


(23) 


Taking the natural logarithm of Eq. (23), we then have, by analogy with Eq. (10), 


/- Hl 

2 VP r > . PT 

— — L _ In 77 4- — 

Nq < A^o 

^0 


(24) 


Since, as previously noted, the term PT /No can be absorbed into the decision threshold, then an equiv- 
alent test to Eq. (24) is 


2 sfp 
No 


L 





Ho 


(25) 


which, is identical to Eq. (11)! Thus, we conclude that in this particular circumstance, the MLR test 
(pseudocoherent receiver ) and the ALR test (noncoherent receiver) are identical Hence, the performance 
of the pseudocoherent receiver is also described by Figs. 2 and 3. It is to be emphasized again that 
the equivalence found here between ALR and MLR receivers is not typical and applies only in this very 
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special case of the detection of a signal with known frequency and unknown phase. More often than 
not, the receiver derived from the MLR approach will have an inferior performance to the optimum one 
derived from the ALR approach. 


IV. A More Precise Formulation of the Problem 

In reality, the subcarriers that are transmitted to indicate the status of the S /C are continuous square 
waves that biphase modulate the carrier. Thus, denoting the carrier radian frequency and phase by u c 
and 9 C (previously called 9 ), respectively, and the square-wave subcarrier radian frequency and phase by 
u sc and 0 SC , respectively, then the received signal analogous to Eq. (la) is given by 

r(t) = s(t , 9 C) 9 SC ) + n{t) = \/2Psin (u c t + 9 C -b ^ Sq ( u sc t -b 0 SC )) + n(t) 


= a/ 2 P Sq (u > sc t + 0 SC ) COS (uf c t + 0c) + n W 


(26) 


Assuming that the harmonics with frequencies other than the sum and difference of u c and u sc are filtered 
out, then in so far as detection is concerned, we may consider the received signal to be 3 

r(t) = s (t,9 c ,9 sc )+n(t) = \[P {cos[(cl> c + u; sc )t + (0 C + 0 SC )\ + cos[(u> c - u sc )t + (6 C - 0 3 c)]}+™(<) (27) 


i.e., the problem is to detect the presence or absence of two tones in an AWGN background where both 
uj c and Use are assumed to be known but both 9 C and 9 SC are assumed to be completely unknown. For 
convenience of notation, we shall rewrite Eq. (27) as 


r(t) = s(£,0+,0_) + n(t) — y/P {cos[u+t 4- 0+] + cos [u-t + 0_]} -b n(t) 


(28) 


where 


U± — uJ c i U sc 

9± = 6 C ± 6 sc 


(29) 


At first glance, it might appear that, because the phases 9 C and 9 SC appear in the two signal tones as 
their sum and difference, the detection of these tones cannot be performed independently. Interestingly 
enough, 9+ = 9 C + 9 SC and 0_ = 9 C — 9 SC when reduced modulo 27T can be shown to be independent 
uniformly distributed random variables (see the Appendix). Thus, as we shall see shortly, the detection 
of two distinct sinusoidal tones with independent random phases in an AWGN background can be treated 
by a likelihood ratio approach analogous to that discussed in the previous section for a single tone in the 
same background. 


* In reality, the >/F amplitude factor in Eq. (27) should be (2\/2/7r)\/P = 0.9003\/P to account for the amplitude of the 
first harmonic in the square-wave subcarrier. For simplicity, we shall ignore this minor difference. 
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A. The ALR Test 

As discussed in Section II, the optimum decision rule is, in general, obtained by applying the average- 
likelihood approach, which in this case means averaging the conditional likelihood function over the two 
random phases 8 + and 0-. In particular, the conditional pdf of the received signal under hypothesis H\ 
is analogous to Eq. (5): 


P{r{t\H,) = /:/; p{r(t)\H u 8 + ,e-)po + ,e_{8 + ,8-)d8 + dd- = J J p(r{t)\H x , 8+, 8-)d8 + d8. 

(30) 


and, hence, the ALR becomes 


A(r{t)) 


2 ^) I-* J*„P{r{t)\H u d + ,8~)d8 + dd_ 


p(r(t)\H 0 ) 


= exp 


r pt] ( i \ 2 r (2 Vp r 

{ Nq J \27tJ /.-( No i 


( 2 Jp r T 1 

r(t) cos (uj-t 4- 9- )dt > dO. 


/ 7T f n / p /*‘J 

exp | ^ — J r(t) cos(uj+t + 0 + )d£ | 




Defining the I and Q correlations for the sum and difference frequencies by 


(31) 


Lc± 


l cos uj±tdt 


= [ r(t)y/2 c 

Jo 

rT 

L s ± = / r(£)\/2sinu;±td£ 


(32) 


then, the likelihood function of Eq. (30) can be rewritten as 




(33) 


where, analogous to Eq. (9), the envelopes corresponding to the upper and lower subcarrier tones are 
given by 


L ± “ + L 2 S± 


(34) 


Alternately, in terms of the log-likelihood function, we arrive at the decision rule 
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(35) 


In J 0 



+ In /o 



H i 


> 

< 


In rj + 


PT 

Ao 




Note that now, despite the fact that In Iq(x) is a monotonic function of its argument, x, we cannot directly 
simplify Eq. (35) to a form analogous to Eq. (11). Rather, to get such a form, one must approximate the 
ln/ 0 (x) function by its series and asymptotic forms for small and large arguments, namely, 


lnl 0 (x) = 


x 

— , small x 
4 

|x|, large x 


(36) 


Thus, for example, if we invoke the small argument approximation of the In Iq(x) function in Eq. (35), 
we get the decision rule (optimum for small SNR) 


Hi 

Ll + L\ 4 L 2 J 7 (37) 

Ho 

where 7 is again a normalized threshold [not necessarily equal to the one defined in Eq. (12)]. The 
decision rule in Eq. (37) suggests the ALR structure illustrated in Fig. 5, which is analogous to that given 
in Fig. 1. For the large argument approximation of the \nIo(x) function, the implementation of Fig. 5 
would require square root devices in each arm entering the final summer prior to the decision device. 

B. The MLR Test 

Let us now again compare the noncoherent two-tone detector derived from ALR considerations and 
specified by the decision rule of Eq. (35) to a pseudocoherent detector that can be derived from MLR 
considerations. In particular, consider the joint ML estimates of defined as 


@ m l+ 1 $ m L— = max 


p(r(t)\H u 6 + ,6_) 


*+.*- p(r(t)\H 0 ) 

which, because of the independence of 6 + and 0_, is determined as 


(38) 


$ M l+ ■> @ m L— = max 
9 + ,0 


f 2 \fp r T ) { 2 s/p f T 

x exp < ^ J r(t) cos (u-t + 6_)dt > exp < — - J r(t) cos(a j+t + 0 4 


)dt 


= »iax |^^L_cos(0_ +a_)|exp|^^L + cos(0 + +a+)| 
where a± are defined in terms of L±, analogous to Eq. (9). The solution to Eq. (39) is 


(39) 


®ML± = 


(40) 
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Fig. 5. Average-likelihood (noncoherent) detector for detection of a pair of independent 
sinusoidal tones with known frequencies and unknown phases in AWGN. 


which, when substituted in Eq. (39), gives 

A(r(())ae x p{-^}exp|^L + |exp|^L.|=exp{-^|exp|^£j (41) 

Taking the natural logarithm of Eq. (40), we get the decision rule 

^(£ + + L_)> M+f (42) 

Ho 

Comparing Eq. (42) with Eq. (35), we observe that, in the two-tone case, the MLR test (which would 
lead to a pseudocoherent form of detector analogous to Fig. 4) is not the same as the ALR test. However, 
using the large argument approximation of the lnio(x) function as given by Eq. (36), we see that the 
ALR and MLR tests once again become equivalent. In summary then, we observe that, for detection of 
a single tone in AWGN, the ALR (noncoherent) test and MLR (pseudocoherent) test are equivalent for 
all SNRs, whereas for the detection of a pair of equal power tones in AWGN, the ALR and MLR tests 
are equivalent only at sufficiently large SNR. 

C. Performance (Receiver Operating Characteristic) 

The performance of the low SNR receiver in Fig. 5 is, as before, described in terms of its false alarm 
probability ( Pp ) and its probability of detection (Pd)- These probabilities are readily computed from 
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knowledge of the first and second moments of the Gaussian random variables L c ± and L s ± [see Eq. (31)] 
under the two hypotheses, namely, 


Ho : E{L C ±} = E{L S ±} = 0 


var { L c ± } = var {L s± } = 


N 0 T 


H\ : E{L c± \0) = \ —T cos 0± 


E{L s± \e} = \/ -^T sin 9± 


var { L c ± } = var {T s ±} 


N 0 T 


(43) 


To compute P F , we observe as before that, under hypothesis Ho, L 2 is a central chi-squared random 
variable (now with two more degrees of freedom). Thus, 


P F = Pr {Hi\Hq} = Pr{L 2 > 7 ^ 0 } = T ^rexp(-r)dr = (l + exp 


(44) 


Similarly, we observe that, under hypothesis H\,L 2 is a noncentral chi-squared random variable (now 
with two more degrees of freedom). Thus, 


P D =Pv{H l \H i }=Pv{L 2 > 1 \H l }= I 


y/27/ NuT \ d 


R ( w ) exp 


R 2 + d 2 \ ^ ( Rd j dR = / 2 7 


N n T 


(45) 


where d 2 is the detection SNR defined as before [see Eq. (16)] and Q]^(a,j3) is the generalized Marcum 
Q-function defined by 


/ z \ ^ ~~ 1 - 4 - Ck^ \ 

Qm(<*,(3) = J z J exp y ) I M -\{otz)dz (46) 

Note that Qm{oi,P) can be obtained from Q(a,f3) = Q\(a,(3) by the relation [2, Appendix 5A] 

Q M (a,P) = Q(a,P) + ex p (°-±^ £ Ij(a0) 

Unfortunately, the normalized detection threshold cannot be explicitly eliminated in Eqs. (39) and (40) 
to give a closed-form expression for the receiver operating characteristic (ROC) analogous to Eq. (18). 
However, for any range of interest, the ROC can be determined numerically. Such numerical results are 
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superimposed on the single-tone detection in Figs. 2, 3(a), and 3(b). We observe that the performance 
penalty associated with using a pair of subcarrier tones each with half the total power relative to the 
full-power single carrier tone case is quite small, e.g., on the order of 0.4 dB or less for P F = 1CT 2 and 
on the order of 0.3 dB or less for P F = 10“ 4 . The degradation associated with the true optimum ALR 
scheme as described by the decision rule of Eq. (35) would be even smaller. Thus, the performance curves 
of the true optimum ALR scheme for two tones would lie between the solid and dashed curves in Figs. 2, 
3(a), and 3(b) since indeed these performance results cannot beat those corresponding to the single-tone 
case. Because of the small degradations involved, we choose not to simulate the true optimum case. 


Part 2. Unknown Frequency and Unknown Phase 
V. The Average-Likelihood Ratio Test 

A. Derivation of the Optimum Decision Rule and Associated Receiver Structure 

Consider the transmission of a fixed (known) amplitude sinusoid with unknown frequency and unknown 
phase over an AWGN channel. Analogous to Eq. (1), the received signal can be modeled over the interval 
of observation 0 < t < T as corresponding to either of two hypotheses, namely, 

7'(t) = s(t , 6) 4- n(t) = V2Pcos(ujt + 6) 4- n(t) (47a) 

when indeed the signal was sent (hypothesis Hi) or 


r(t) = n(t) 


(47b) 


when the signal wets not sent (hypothesis Hq). In addition to the previously defined parameters, in 
Eq. (47a), / = lu/2i r denotes the unknown carrier frequency assumed to be uniformly distributed in the 
interval (f c — B/2, f c + B / 2), where as before f c denotes the nominal carrier frequency (i.e., in the absence 
of Doppler). When the signal has two unknown parameters, e.g., the phase 0 and frequency /, then to 
compute the numerator of Eq. (3), we must first condition the pdf p(r(t)\H\) on both of the unknown 
parameters and then average over them, i.e., 

rfrA-B/2 rTT 

p(r(t) \Hi)=l / p(r(t)\H l) ej)p e (9)p f (f)dedf (48) 

Jfa-Bf 2 J-ir 

where Pe(&),Pf(f) respectively denote the pdf’s of the unknown parameters 6 and /. In our situation, 
the phase and frequency are assumed to be completely unknown, and thus pe(6) and p/(f) are uniform 
distributions. Hence, combining Eqs. (3) and (48), the average-likelihood ratio (ALR) becomes 


A(r(t)) 


1 

2nB 



p(r(t)\H u e,f)d6df 


p(r(t)\H 0 ) 


exp 


exp 


f_PT| J_ rf * +B ? 2 

I Wo / 27 tB Jf' -B/Z 

\ PT] 1 r^ +B / 2 

iNo/S J fc _ B/2 0 



4v/p 

No 


L(f)J df 



(49) 


83 



where 


L(f) 4 y/L\{f) + L\(f) 


with 


rT 

L c (/) = / r(t)\/2cos27Tftdt 

Jo 

L s (f) = J r(t)y/2 sin ft dt J 


(50) 


(51) 


It should be noted that L(/) is nothing more than the magnitude of the complex Fourier transform (FT) 
of r(t) in the interval 0 < t < T. If r(t) is band limited to W Hz, then for large WT , the real and 
imaginary components of this complex FT, namely, L c (f), L s {f) can be approximated by the discrete 
Fourier transforms (DFTs) 


- 2 WT 

LM = " ( 2 #) C “ 

n= 1 
2 WT 

i.(/) = Y, r ( 2 ^) siu ( 2 i /5w) 


(52) 


Comparing A {r(t)) to a threshold produces (after suitable normalization) the ALR decision rule 



(53) 


Since Eq. (53) is overly demanding to implement, one discretizes the frequency uncertainty interval into 
G = B/T~~ l - BT subintervals to each of which is associated a candidate frequency /*; i — 0, 1 , 2, * • • , G-\ 
located at its center. As such, the integration over the continuous uncertainty region in Eq. (53) is 
approximated by a discrete (Riemann) sum and, hence, the approximate decision rule becomes 


G - 1 


s> 



Hi 

> 

< ^ 
Ho 


(54) 


which has the implementation representation of Fig. 6. It is important to understand that spacing the 
frequencies 1 ~ 0, 1, 2, • • • , G — 1 by 1/T guarantees independence of the noise components that appear 
at the output of each spectral estimate channel. However, orthogonality of the signal components of these 
same outputs depends on the true value of the received frequency relative to the discretized frequencies 
f , \ i = 0, 1, 2, ■ • ■ , G— 1 assumed for implementation of the receiver. That is, if the true received frequency 
happens to fall 011 one of the fi s, then a signal component will appear only in the corresponding spectral 
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Fig. 6. Average-likelihood (noncoherent) detector for detection of a single sinusoidal tone with unknown frequency and 

unknown phase in AWGN. 


estimate channel, i.e., all other channels will be noise only. On the other hand, if the true received 
frequency falls somewhere between two of the ff s, then we have a loss of orthogonality in that a spillover 
of signal energy occurs in the neighboring spectral estimates. The worst-case spillover would occur when 
the true received frequency is midway between two of the ff s. 

We conclude this section by noting that a decision metric similar to Eq. (54) arises in the study of 
FH or DS/low probability of intercept (LPI) optimum ALR (noncoherent) detection [3-5], where in the 
FH case, /*; i — 0, 1, 2, ■ * • , <3 — 1 corresponds to the G possible frequencies that the transmitted signal 
can hop to and the detection is based on observation of a single hop of duration Th = T , and in the DS 
case G is the number of possible code sequences that can occur in the observation interval. Many of the 
results obtained from these works are directly applicable to the problem at hand. 

B. Performance 

It is tempting for large values of G (as is typically the case) to apply a central limit theorem argument 
to the left side of Eq. (11), i.e., approximate it as a Gaussian random variable in so far as computing the 
receiver operating characteristic associated with this decision rule [4]. Unfortunately, it was shown in [5] 
that following such an approach is very poor when compared with results obtained from simulation or 
numerical methods applied to the true decision rule of Eq. (11), even for values of G as large as 1000 or 
10,000. In fact, it is stated in [5] that G on the order of “ten thousands is not guaranteed to be large 
enough to validate the Gaussian approximation.” Thus, to obtain the true receiver performance, we too 
must resort to simulation and/or numerical methods, such as those suggested by Requicha [7], wherein 
the characteristic function and fast Fourier transforms (FFTs) are used to compute approximate values 
of the distribution function associated with the left-hand side of Eq. (11). More about this later on. 
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VI. The Maximum-Likelihood Ratio Test 

A. Derivation of the Optimum Decision Rule and Associated Receiver Structure 

Although the exact evaluation of the numerator of the likelihood ratio in Eq. (2), i.e., p(r(t)\H\) is 
obtained from the law of conditional probability as described by Eq. (5), namely, conditioning on the 
unknown parameters and averaging over their distribution, it is also possible to approximate this numer- 
ator by first finding the ML estimates of the unknown parameters and then substituting these values into 
the conditional probability p(r(t)\Hi, 0, /). That is, we approximate p(r(t)\Hi) by p(r(t)\Hi,Q ml, f ml), 
in which case the likelihood ratio test (now referred to as the ra^mmura-likelihood ratio (MLR) test) 
becomes 


A (r(t)) * 


P (r{t)\HuOb xl, /ml) 
p(r(t)\H 0 ) 


H i 
> 
< 
H 0 


where 


Oml, f ml = max 
vj 


p(r(t)\H 0 ) 


(55) 


(56) 


The maximization over 6 required in Eq. (56) can be performed identically to that in Section III [see 
Eq. (23)]: 


max 

9 


P(r(t)\H 

P(r(t)\H 0 ) 


= exp 




(57) 


where L(f) is as defined in Eq. (50). Thus, the optimum maximum a posteriori (MAP) decision rule 
becomes 


max exp 
/ 



H i 

> 

< 71 
H o 


(58) 


Since the exponential is a monotonic function of its argument, we have the equivalent decision rule 4 


H i 

max L(f) > y/j (59) 

Ho 

which results in a spectral m.axim,um form of receiver. Again, because of the excessive demand placed on 
the implementation by the need to evaluate Eq. (59) over a continuum of frequencies, we again quantize 
the frequency uncertainty interval into G = BT subintervals, each with an associated candidate frequency 


4 We define the normalized threshold equal to ^ to be consistent with the notation used in Part 1. In this way, when G 
is equated to unity, then our results obtained here will reduce to those given in Part 1 for the MLR decision rule. 
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fu i — 0, 1, 2, ■ ■ , G - 1 located at its center. Thus, the frequency continuous decision rule of Eq. (59) can 
be approximated by the decision rule 


H i 

maxL(/i) < V7 (60) 

Ho 

which suggests the receiver of Fig. 7. Here again, as with Eq. (54), the orthogonality of the spectral 
estimates is not guaranteed unless the frequency of the received signal falls on one of the /Vs. Also, since 
£(/;); * = 0, 1 , 2, • • • , G represents a uniform sampling of !(/), then in view of Eq. (52), we can implement 
Fig. 7 with FFT techniques. 



Fig. 7. Maximum-likelihood detector for detection of a single sinusoidal tone with unknown frequency and 

unknown phase in AWGN. 


B. Performance 

The performance of the MLR decision rule of Eq. (60) can be obtained analytically since the pdf of 
G independent random variables can be explicitly written in terms of the pdf’s of individual random 
variables, which in turn are obtained from the results in Part 1. The procedure is as follows. 

1. Best-Case Performance. Consider first the optimistic (best) case, where the actual received 
carrier frequency is indeed equal to one of the G frequencies, say fk used to approximately implement 
the optimum decision rule as per the discussion following Eq. (59). Under H\,G — 1 of the L(/i)’s are 
Rayleigh distributed with pdf [see Eq. (14)] 

= £>o ( 61 ) 
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and the single L(f l ) that is associated with the received signal carrier frequency, namely L(/fc), is Rician 
distributed with pdf [see Eq. (15)] 


PW „m = ^exp(-^)/„(^) ; L>0, s-pi* («> 

Let Pp denote the per frequency channel false alarm probability, i.e., 

which is independent of fi. Then, the overall false alarm probability, Pp, is given by 
P F = Pr | max L(/j) > V7l-^o} 

- 1- Pr {L(fo) < s/l-, L(fi) < y/i, • • • , L(fc-i) < y/l\Ho} 

G - 1 G - 1 

= 1 - n Pr W‘) ^ Vl\ H o} = 1 - n 0 - Pr W*) > yfr\H 0 }) 

1=0 1=0 

.l-(l-p?) G =l-(l-«p (64) 

Since, under Hi, G — 1 of the spectral estimates (i.e., the ones containing noise only) have the same pdf, 
namely Eq. (61), as under Hq , and one spectral estimate has the Rician pdf of Eq. (62), then the overall 
probability of detection, P d , is determined from 

G - 1 

P D = Pr | max L(f l ) > = 1- (1 - Pr{L(/;) > y/j\H],}) = 1-(1 - P F ) G 1 (1 - Pp) ( 65 ) 

1=0 

where Pb corresponds to the detection probability of the 
i.e., 


P' D = Pr {£,(/,)> s /j\H l } = Q^d. 
Substituting Eqs. (63) and (65) into Eq. (66) gives 

p D = (-^)y 

or, equivalently, the overall probability of miss, P F j, is 


single-frequency channel containing the signal, 


27 \ , 2 2PT 

NoTj ’ N 0 


( 66 ) 




( 67 ) 
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P " 4 l - Pc= (l-e x p(-A^)) 


(68) 



Note that, for G = 1, Eqs. (64) and (67) reduce, respectively, to Eqs. (14) and (15). 

The ROC can be determined by eliminating the normalized threshold between Eqs. (64) and (67), in 
which case one obtains 


Pd = 1 - (1 - Pf) {G ~ 1)/G (l-Q (d, y/-2]n(l-(l-P F )'/G)y) 


(69) 


2. Worst-Case Performance. The worst-case performance occurs when the actual received carrier 
frequency is indeed midway between two of the G frequencies used to approximately implement the 
optimum decision rule as per the discussion following Eq. (59). Under Hq, the false alarm performance is 
still described by Eq. (64). However, under H i, all G spectral estimates are now Rician distributed with 
pdf’s of the form in Eq. (62), namely, 


PLUi){L ) - 


2 

No T 


Lexp 



L > 0 


where the /Vs are determined as follows. Since [see Eq. (15)] 


(70) 


0i = (E {L e (fi)\9, f}f + (E {L a (fi)\e, f }f 


(71) 


then, assuming that the actual received carrier frequency, /, is situated midway between fa and fa+ 1 , 
which are separated by 1/T, i.e., / = f k -f 1/2T, Eq. (70) is evaluated as (for simplicity, we ignore the 
edge effects at the ends of the frequency uncertainty band) 



i — fc, k + 1 
i ^ k,k 4-1 


= pr 2 r t 

Finally then, analogous to Eq. (70), the detection probability would be given by 


(72) 


p D = i- n (i - Q ( r i 

i=0 ' ^ 



(73) 


which, in general, depends on /^, i.e., the location of / within the uncertainty band. 

It has been suggested in [5] that the two nearest spectral estimates (envelopes) to the frequency 
location of the received signal dominate the performance, i.e., the spillover effect of signal in the other 
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frequency slots can be ignored to a first-order approximation. When this is done, then, under H\, two of 
the spectral estimates are identically Rician distributed and the remaining G - 2 are identically Rayleigh 
distributed. In this case, Eq. (73) is replaced by an expression somewhat like Eq. (67), namely, 


Pd = 1 



(74) 


which is now independent of the frequency location of the signal. Combining Eqs. (64) and (74), the 
ROC is approximately given by 


P D = 1 - (1 - P F ) iG 


— 2 )/G 


1 - Q ((f) d * >/ _ 21n(l - (1 -Pf) 1/g 


(75) 


VII. A More Precise Formulation 

As discussed in Section IV of Part 1, the true transmitted signal corresponds to a sinusoidal carrier 
phase modulated by a square- wave subcarrier of radian frequency uj sc . At the receiver, the harmonics 
with frequencies other than the sum and difference of tu sc and lj c are filtered out, which means that in so 
far as detection is concerned, the received signal in the absence of frequency uncertainty can be modeled 
as 


r(t) = s(t,6 c , 6 sc ) + n(t) = y/P { cos [(cj c +c o sc )t -f (6 C + <9 SC )] + cos [(w c - co sc )t + ( 6 C - <9 SC )]} + n(t) 

(76) 


In the presence of frequency uncertainty due, for example, to Doppler shift, both the upper and lower 
frequency tones in Eq. (76) will be shifted from their nominal values with the higher-frequency tone 
experiencing a larger shift than that corresponding to the lower-frequency tone. If, however, the subcarrier 
frequency is much smaller than the carrier frequency, i.e., lj sc << u c , as is the case of interest, then for all 
practical purposes, one can associate the frequency uncertainty with the carrier as discussed in Section V. A 
and assume to a first-order approximation that both upper and lower frequency tones experience the same 
frequency shift. Stated another way, we can assume that, in so far as detection is concerned, we observe 
a pair of tones whose frequencies are unknown (but by the same amount), each in a band D Hz centered 
around its nominal value. Furthermore, the uncertainty band is assumed to be very narrow with respect 
to the subcarrier frequency, i.e., B « f sc . 

A. The ALR Test 

Analogous to what was done in Part 1, the conditional pdf of the received signal under hypothesis H\ 
is given by 


p(r(t) m 




p(r(t)\H u 6 + , 0_, / - / sc , / + f sc )d6+dd-df 


(77) 


whereupon the ALR becomes 
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A(,«» = ex P {^} 5 «t 


(78) 


In Eqs. (77) and (78), the spectral envelopes at the lower and upper tones are defined by 


L±(f) = 


( 79 ) 


together with 


L c ±{f) = [ r(t)V 2 cos[ 2 Tr(f ± f sc )t]dt 

Jo 

L s ±(f) = J r (t)V 2 sin [ 2 ir(f ± f sc )t]dt j 


Discretizing the integration interval results in the approximate decision rule 




X H ' 

L-(fi ) J 
' Ho 


(80) 


(81) 


where the spectral envelopes required in Eq. (81) are defined analogously to Eqs. (79) and (80), with the 
continuous random variable / replaced by the discrete random variable /*; i = 0, 1, • • • , G — 1. As was 
the case for the single-tone result in Section V.A, the performance (ROC) of the decision rule in Eq. (81) 
cannot be obtained analytically. 

B. The MLR Test 

Without going into great detail, it is straightforward to show (using the results of Section IV. B) that 
the MLR test analogous to Eq. (58) becomes 


T exp exp (V M/) ) exp (‘ W LM) ) < ” 

Ho 

or, equivalently, 

Hi 

ma x{L-(f) + L+(f)) > ^7 
Ho 

which has the discretized version 

Hi 

max(L_(/i) + £+(/)) < a/7 

l i 

Ho 


(82) 


(83) 


(84) 
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Unfortunately, the performance of the receiver that implements the decision rule of Eq. (84) also cannot 
be obtained analytically. 

VIII. Numerical Results 

Since the performance of none of the ALR optimum decision rules can be evaluated analytically and 
since the same is true for some of the MLR decision rules, a computer simulation of these metrics has 
been developed to numerically evaluate such performance. The results of such simulations are described 
as follows. Figure 8 is a sample illustration of the ROC for the case of a single tone with unknown phase 
and frequency (as described in Section V) and a detection SNR d 2 — 2PT/Nq = 6 dB. Both ALR and 
MLR cases are illustrated, corresponding, respectively, to the decision rules of Eqs. (54) and (60). Also, 
both the best- and worst-case input frequency scenarios are considered, corresponding, respectively, to the 
cases where the actual input frequency is indeed equal to one of the G frequencies used to approximately 
implement the decision rule and the case where the actual input frequency falls midway between any 
two of these G frequencies. Clearly, the actual system performance corresponding to an input frequency 
arbitrarily chosen in the uncertainty band will lie between these two performance bounds. We observe 
from the results in Fig. 8 that the difference between best- and worst-case performance is relatively small, 
as well as is the difference between the ALR (optimum) and MLR (suboptimum) decision rules. There 
is a significant difference, however, between the performance for G = 10 and G = 100, indicating the 
sensitivity of the performance degradation to a factor of 10 increase in frequency uncertainty. Also, 
comparing Fig. 8 with the analogous curve in Fig. 2, corresponding to the case of unknown phase but 
known frequency, we again see a rather significant degradation in performance when the frequency is 
unknown even by only a factor of 10 relative to the observation bandwidth (reciprocal of the observation 
time, T), i.e., G = 10. 



Fig. 8. ROC: frequency and phase unknown (single-tone) 
simulation results. 


As verification of the MLR simulation results, we present in Fig. 9 the analogous analytical results 
obtained from Eqs. (69) and (75). Recall that in arriving at Eq. (75) the assumption was made that 
t he energy spillover effect of the signal into the other frequency slots is dominated by the two adjacent 
ones. Thus, ignoring edge effects, it was not necessary to average over all possible worst-case (mid- 
way) input frequency positions. In the computer simulation, this assumption was not invoked, as the 
input frequency was allowed to occur midway between any two adjacent frequencies. Despite this analysis 
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approximation, however, comparison of the results in Figs. 8 and 9 reveals excellent agreement between 
analysis and simulation, i.e., the assumption of only adjacent signal energy spillover used to arrive at 
Eq. (75) has been justified. Also indicated in Fig. 8 is the analytical result corresponding to known phase 
and frequency (recall that this result is the same for both MLR and ALR) that allows a more direct 
assessment of the performance degradation due to lack of perfect frequency knowledge. 

Since the curves in Fig. 8 are drawn for a fixed value of detection SNR d 2 = 2PT/N 0) then assuming 
that P/N 0 is specified, this implies that the observation interval, T, is also held constant. Thus, changing 
the value of G = BT from 10 to 100 directly translates into a change by a factor of 10 in the frequency 
uncertainty region B , which accounts for the observed degradation in performance. Another interpretation 
of the numerical data can be obtained by again holding P/No fixed but observing the effect on system 
performance of increasing T for a fixed frequency uncertainty region B. This necessitates plotting the 
ROC with both d 2 and G increasing linearly with T. Such a plot for the ALR decision rule with best-case 
input frequency is illustrated in Fig. 10, where the ROC is plotted for values of G ~ 10, 20, 40, and 80 (T 
increasing by a factor of 2) and corresponding values d 2 = 6, 9, 12, 15 dB. To directly see the dependence 
of MLR system performance on detection SNR, Fig. 11 illustrates the behavior of detection probability, 
Pd, versus detection SNR, d 2 , for a fixed false alarm probability, Pp = 10 -2 , and values of G — 10 and 
100. These curves are obtained from numerical evaluation of the analytical results in Section VI. Since 
along any curve G is held fixed, one can interpret these results as keeping the frequency uncertainty band, 
R, and observation time, T, fixed and observing the change in performance as P/No is varied. 

The penalty associated with detecting a pair of subcarrier tones (each at half the total transmitted 
power) as opposed to a single carrier tone (at full transmitted power) is illustrated by the numerical 
results in Fig. 12. Here we plot the ROC for both the single- and double-tone cases for the ALR decision 
rule with best-case input frequency and a detection SNR equal to 6 dB. The results for the single-tone 
case are taken directly from Fig. 8. We observe a significant performance penalty associated with using a 
double- tone detection scheme. Figure 13 illustrates for the double- tone detection scheme results analogous 
to Fig. 10 for the single-tone detection scheme. Here again, by comparing the two figures, we observe a 
significant penalty associated with using a pair of equal half-power subcarrier tones rather than a single 
tone at full power. 



Fig. 10. ROC simulation results: frequency and phase unknown 
(single tone), ALR , best-case input frequency. 
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Fig. 13. ROC simulation results: frequency and phase unknown 
(double tone), ALR, best-case input frequency. 
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Appendix 

On the Independence of the Sum of Difference of Two Uniformly 
Distributed Random Variables Modulo 2n 


Consider two independent random phases 9 A and 0 B that are each uniformly distributed in the semi- 
closed interval [— 7T, 7 t). Define the sum and difference of these two random variables by 


o' + =o A +e B ' 
o'_ =e A -e B , 


and the modulo 2n versions of these random variables by 


9+ - {9'+) mod 27r - (9a + 9b ), nod 2x 
«- = (*-)„* 2, “ («A - M mo d 2. 


(A-l) 


(A-2) 


The probability density functions (pdf’s) of 0' + and 6'_ are triangular in the semiclosed interval [— 27 t, 27 r), 
i.e., they are the convolutions of two uniform pdf’s, whereas the pdf’s of their modulo 27 t reduced versions, 
0+ and 0_, are once again uniformly distributed in [-7r,7r) (see Fig. A-l). We would now like to show 
that 0+ and 6 _ are indeed independent random variables. To do this, we shall show that the conditional 
pdf P0_(9-\O + ) satisfies p0_(0_|0+) = pe_(9 _), i.e., it is a uniform distribution in [— 7r,7r). Similarly, it 
can be shown that p0 + (O+\O _) — p# + (0 + ). 

Let 9+ be any positive value in its region of definition, i.e., 0 < 0+ < 7r. Then, 0 A and 0 B are related 
as follows: 


A _ f -0A 4- 0 + - 2i r, -7T < 9 a < -7T + 

B “ \ ~0 A + 0 + , -7T + 0 + < 9 A < 7T 

From Eq. (A-l), we find that 

_ j 29 A - 0+ 4* 2?r, -7t < 9 a < -7r + 0+ 

^ 29 a - 0+, -7T + 9 + < 9 a < 7i 


(A-3) 


(A-4) 


Thus, from Eq. (A-4) and the fact that 9 A is uniform in the interval [ — 7r, 7 t), the conditional pdf p e >_ (0'_ |0-j-) 
appears as in Fig. A-2(a). Reducing 0+ modulo 2tt produces the conditional pdf p0_ (0_ |0 + ) as illustrated 
in Fig. A-2(b), i.e., a uniform distribution in the interval [-7r,7r) Q.E.D. 
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